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Abstract 

The study of DaUtz plots of heavy-meson decays to multi-hadron final states has received intensified interest by the 
possibility to gain access to precision investigations of CP violation. A thorough understanding of the hadronic final- 
state interactions is a prerequisite to achieve a highly sensitive, model-independent study of such Dalitz plots. We 
illustrate some of the theoretical tools, predominantly taken from dispersion theory, available for these and related 
purposes, and discuss the low-energy decays o), (p 3n in some more detail. 
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1. Dalitz plots and CP violatioifl 

A precise study of final-state interactions is increas- 
ingly becoming of high importance for our understand- 
ing of the most diverse aspects of particle decays in- 
volving hadrons. Final-state interactions can be of sig- 
nificance for various reasons: if they are strong, they 
can significantly enhance decay probabilities; they can 
significantly shape the decay probabilities, most promi- 
nently through the occurrence of resonances; besides 
resonances, also new and non-trivial analytic structures 
can occur, such as threshold or cusp effects (for the 
prominent role cusp effects have played recently in 
studying pion-pion interactions, see Ref. |2] and refer- 
ences therein); and finally, they introduce strong phases 
or imaginary parts, the existence of which is a prereq- 
uisite for the extraction of CP-violating phases in weak 
decays (see e.g. Ref. |3]). Dalitz plot studies of weak 
three-body decays of mesons with open heavy flavor 
(both D and B) are expected to acquire a key role in 
future precision investigations of CP violation, due to 
their much richer kinematic freedom compared to the 
(effective) two-body final states predominantly used to 
study CP violation at the B factories. In many cases. 
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the branching fractions are significantly larger; further- 
more, the resonance-rich environment of multi-meson 
final states may help to enlarge small CP phases in parts 
of the Dalitz plot, and differential observables may al- 
low to obtain information on the operator structure that 
drives CP violation beyond the Standard Model, once 
it is observed. Since the results from the B factories 
have shown that the Cabibbo-Kobayashi-Maskawa the- 
ory 10] represents at least the dominant source of CP vi- 
olation, our long-term goal will be to find other sources 
of CP violation that contribute additional, smaller ef- 
fects. For this purpose, clearly extremely accurate mea- 
surements and means of theoretical interpretation are 
required. Strong evidence for CP violation in three- 
body final states has already been reported for B^ 
K^n^Ti^ (3], e.g. with a 3.7cr signal in the effective Kp 
channel. Only negative results exist for D decays so 
far 10] ; however with the first preliminary evidence for 
direct CP violation in D° K+R- vs. D° ^ tt+tt" 
reported [7J, cross checks of any mechanism proposed 
for its explanation in three-body decays would be of 
paramount importance. 



There are different possibilities how to analyze CP 
violation in Dalitz plots. One suggestion is a strictly 
model-independent extraction from the data directly 
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1^, e.g. using the significance variable |@] defined as 



^JN(i) + Nil) 



(1) 



where A^(/) and N(i) denote the event numbers of CP- 
conjugate decay modes in a specific Dalitz plot bin 
CP violation can then be identified in a deviation from 
a purely Gaussian distribution in the significance plots. 
The significance method allows to study local asymme- 
tries and requires no theoretical input at all. 

An alternative approach is to make use of as much a 
priori theoretical information on the strong amplitudes 
as possible. To see why this may be advantageous, con- 
sider the following toy modelj^ let the event number dis- 
tributions be given by a (Breit-Wigner) resonance sig- 
nal (of mass Mres and width r,es, with strength /3) on 
a certain background (oc a), with a CP-violating phase 
6cp, according to 



N, N = a-Hj6Re 



!,±1*CP 



S - Ml, + /Mi-esFr, 
sin 6cp X 2jeMresrres 
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(2) 



As the pole position of the resonance in question (given 
by Mies and Ties) is universal and may, for the sake of 
the argument, be accurately determined by other, inde- 
pendent processes, it is obvious that the functional form 
of Eq. (|2]i is strongly restricted. Indeed, one can easily 
construct examples where the significance distribution 
shows no deviation from a Gaussian due to lack of data 
statistics [1], while a fit with the strong amplitude (|2]l 
still allows to extract 6cp with some (limited) accuracy 
even in a sparse sample. So while the theoretical as- 
sumptions and prejudices going into such an analysis 
clearly have to be very carefully judged, their benefit in 
terms of vastly increased sensitivity is also obvious. 

In the following, we will therefore briefly sketch 
some of the tools available, mainly based on dispersion 
theory, to analyze the hadronic amplitudes of the (light) 
final-state particles (such as pions and kaons): (light) 
meson scattering, form factors for the description of two 
strongly interacting particles, and the consistent treat- 
ment of three-body final states. 

2. Scattering and form factors 

Analyticity, unitarity, and crossing symmetry pro- 
vide a high degree of constraint for the pion-pion scat- 



disc 



Figure 1 : Graphical representation of the consequence of analyticity 
and unitarity for form factors. 



tering amplitude. They can be exploited using dis- 
persion relations, which can be formulated as a cou- 
pled system of partial-wave equations, the so-called 
Roy equations OlOll . Modern precision analyses of the 
Roy equations have been performed partly mak- 
ing use of constraints from chkal perturbation theory 
on the scattering lengths appearing as subtraction con- 
stants therein fl2], and a similarly rigorous study exists 
also for pion-kaon scattering 1 1,31 . These provide us 
with high-precision parameterizations of the most rele- 
vant scattering amplitudes for light mesons appearing in 
the final states of heavy-meson decays. 

Final-state interactions between only two strongly in- 
teracting particles can be described in terms of form fac- 
tors, which in turn can be linked to the properties of 
scattering amplitudes using analyticity and unitarity. As 
illustrated in Fig.[Tl the unitarity relation for a form fac- 
tor Fj{s) (here: of the pion) of isospin / and angular 
momentum J reads 

Im F'j(s) = F'j(s) X 6{s - AMI) x sin 6'j{s)e-"^'j^'^ , (3) 

from which one immediately deduces Watson's final- 
state theorem lll4ll : the form factor shares the phase 
S'j{s) of the (elastic) scattering amplitude. The solu- 
tion to Eq. (O is obtained in terms of the Omnes func- 
tion [[ll]. 



F'j(s) = P'j{s)D.'j{s) 
- 

^ J41 



AMI ~ 



(4) 
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where P'j{s) is a polynomial. Note that the Omnes func- 
tion is completely given in terms of the phase shift. A 
classic application of such a form factor representation 
is the pion vector form factor Fy{s), making use of the 
nn P-wave phase 6\{s); for a precision parameteriza- 
tion up to about 1 GeV, p - oj mixing and the onset 
of (relatively week) inelasticities have to be accounted 
for il^ ITtIi . Such form factor representations can be 
used for analyses of the e^e' n^n^ data to reduce the 
error of the hadronic contribution to the muon g-2, or to 
check the compatibility of the data with analyticity and 
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Figure 2: Example for the complication of the analytic structure of 
4-point functions through crossed-channel effects, here for the decay 
of a heavy meson (double line) into three pions (single lines). 



unitarity fl?']. Note that the effects of chiral dynamics 
are particularly important for scalar form factors, where 
a parameterization in terms of Breit-Wigner resonances 
can lead to completely wrong phase motions (see e.g. 
Ref. ilST for the context of B — > 3;r decays). 



3. Dispersion relations for three-body decays: 
6>, 3;r 

The application of dispersion relations to three-body 
decays is more complicated than the treatment of form 
factors due to the more involved analytic structure, 
and the possibility of crossed-channel rescattering; see 
Fig.|2]for a depiction of the complication of the unitarity 
relation. One specific (low-energy) example that has re- 
ceived much renewed attention recently lll9ll is 77 ^ 3n, 
due to its importance for the extraction of the light quark 
mass ratios. Here, we concentrate instead on the even 
simpler three-pion decays of the lightest isoscalar vec- 
tor mesons, o), (p ^ ti^ti^tC lEoll . One starts by decom- 
posing the amplitude N\{s, t, u) according to 

M{s, t, u) = e^yapn^'pl^pl-pfla T(s, t, u) , (5) 

where is the polarization vector of the decaying w/^. 



and s = (p,:- + p„+), t = {p„- + p„oy 



iPn* + p^r 



Due to Bose symmetry, only partial waves of odd an- 
gular momentum contribute; neglecting discontinuities 
of F- and higher partial waves, Tis, t, u) can be further 
decomposed as 



r{s, t, u) = T{s) + r(t) + r(u) . 



(6) 



The unitarity relation for Tis), assuming elastic final- 
state interactions, then leads to the following expression 
for the discontinuity of fis): 



disc Tis) ^ 2i{r{s) + f{s)] 

xd{s-4Ml)x sin 6\{s)e 
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(7) 



Equation (|7]i is complicated compared to Eq. (O by in- 

homogeneities 'F{s), which are given by angular aver- 
ages over T according to 



{z"f}is) 
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where A(x, y, T^-x^+y^+z^ - 1{xy + xz + yz), and 
My is the mass of the decaying vector meson. Note 
that the angular integration including the k{s) function is 
non-trivial and generates a complex analytic structure, 
including three-particle cuts due to the fact that uj and (p 
are unstable and decay. For details on how to properly 
deform the angular integration path, see e.g. Ref. 112 ill . 
The analog to the Omnes solution (HJi are then integral 
equations involving the inhomogeneities (compare also 
Ref. |22]) 



ir(^) = Q|(Ja+- r 

\ ^ J AMI S' 



ds' sin6\{s')f{s'y 



(9) 



\0.\{s')\{s' - s) , 

with the subtraction constant a. The number of subtrac- 
tions is chosen such that the dispersion integral is guar- 
anteed to converge. (See Ref. [231 for earlier, related 
formulations.) Again, care has to be taken when per- 
forming the dispersive integral, as 'F{s) shows singular 
behavior at the pseudothreshold s - (My - M„)^. 

Equations (O and (|9]l can be solved iteratively: start- 
ing from an arbitrary input function T{s), we can calcu- 
late the inhomogeneity 'F{s) according to Eq. (jSj, from 
which a new Tis) is obtained from Eq. (|9]l, provided 
we can devise a method to determine the subtraction 
constant; the procedure is stopped once a fixed point 
of the iteration is reached with sufficient accuracy. In 
the example discussed here, see Eq. (|9]l, the subtrac- 
tion constant works as an overall normalization factor of 
the solution; we match it to the partial decay width, but 
note that a normalized Dalitz plot distribution is subse- 
quently a pure prediction. While the result is indepen- 
dent of the starting function, for the case at hand, we 
choose 'F(s) - 0[(s) in order to allow us to quantify 
crossed-channel effects (generated by the iteration) in a 
plausible way. 

Figure [3] shows the result of such an iteration for the 
decay (p — » 37r: it converges fast, with the third iteration 
already all but indistinguishable from the final result. 
The difference to the starting point of the iteration, the 
Omnes function without any crossed-channel rescatter- 
ing, is however very significant. The picture for w — » 37r 
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Figure 3: Successive iteration steps of real (upper panel) and imagi- 
nary (lower panel) part of the amplitude "Fis) for (p 3n. The dashed 
lines denote the physical region of the decay. 



(not stiown here) is qualitatively very similar, with con- 
vergence reached even faster (after two iterations). The 
resulting Dalitz plots for both w — > Stt and — > 37r are 
shown in Fig. ID normalized by the P-wave phase space 
factor (such that what is effectively shown is |;F(5, t, u)\^ 
in arbitrary normalization), using the kinematical vari- 
ables 



2Mv 



y ■ 



{My - M„f - s 
2Mv 



(10) 



The w — > Stt Dalitz plot distribution is relatively 
smooth, rising from the center to its borders, with a 
maximum enhancement in the outer corners of about 
20%. The (/> — > 37r one displays significantly more struc- 
ture, as the p resonance bands lie inside the physical de- 
cay region: the Dalitz plot distribution again rises from 
its center to the three p bands, before falling off steeply 
towards the corners. The profile of the — > 3;7r Dalitz 



plot agrees quaUtatively with experimental results Il24ll : 
a detailed comparison is in progress f2^. 

In order to illustrate the effects of crossed-channel 
rescattering explicitly, we compare the full dispersive 
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Figure 4: Dalitz plots for at —> irr (top) and i 
malized by the P-wave phase space factor. 



37T (bottom), nor- 



solution to the approximation with inhomogeneities 
F{s) in Eq. ^ neglected. By keeping the subtraction 
constant a fixed (assuming, say, it can be determined us- 
ing some other theoretical insight), we find that crossed- 
channel effects enhance the u) —> 37t partial width by 
about 20%, while the ^ —> 3n partial width, in con- 
trast, is reduced by about the same amount. As far as 
the Dalitz plots are concerned, this effect can however 
largely be absorbed in a redefinition of the subtraction 
constant. The remaining effect on the normalized dis- 
tributions is demonstrated in Fig.|5] crossed-channel ef- 
fects increase the Dalitz plot densities in the center, and 
decrease them towards the borders, where the size of the 
effects is significantly larger for the decay of the i^; the 
p resonance bands in the case of —> Stt are left rather 
untouched. 

The method sketched here is currently also applied 
to other light-meson decays such as rj' —> rjnn 11251] . 
Challenges to be faced when extending this formalism 
to heavy-meson decays include the necessity to treat 
systems of integral equations when coupled channels 
within one partial wave cannot be ignored, or inelas- 
ticities are not negligible. In particular when consid- 
ering B-meson decays, elastic unitarity will surely not 
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Figure 5: |!Ffuiil^/|7^^oP for oj ^ 3n (top) and -> 3;r (bottom). 
The subtraction constant is fixed so that it reproduces the decay width 
before and after the iteration. 

be sufficient. It will have to be checked when a purely 
perturbative treatment of crossed-channel effects is fea- 
sible (compare e.g. Ref. |0]), and when higher partial 
waves become important. To investigate these and re- 
lated questions is part of the program of the informal 
Les Nabis network ll27ll . which aims at optimizing fu- 
ture Dalitz plot studies along the lines sketched here — 
with the strong goal to better interpret the mechanism of 
CP violation in nature, yet at the same time learning im- 
portant lessons on nonperturbative strong interactions. 
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